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Geometrical shock dynamics, also called CCW theory, yields approximate equations for 
shock propagation in which only the conditions at the shock appear explicitly; the post- 
shock flow is presumed approximately uniform and enters implicitly via a Riemann in- 
variant. The nonrelativistic theory, formulated by G. B. Whitham and others, matches 
many experimental results surprisingly well. Motivated by astrophysical applications, we 
adapt the theory to ultra-relativistic shocks advancing into an ideal fluid whose pressure 
4—> . is negligible ahead of the shock, but one third of its proper energy density behind the 

shock. Exact results are recovered for some self-similar cylindrical and spherical shocks 
with power-law pre-shock density profiles. Comparison is made with numerical solutions 
of the full hydrodynamic equations. We review relativistic vorticity and circulation. In an 
ultrarelativistic ideal fluid, circulation can be defined so that it changes only at shocks, 
qq ' notwithstanding entropy gradients in smooth parts of the flow. 
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Gamma-ray-burst afterglows have spurred us to look into this problem. These cosmo- 
logically distant events are believed to involve a shock launched by the death of a massive 
star with initial Lorentz factor Tq > 10 2 relative to a pre-shoc k circumstellar wind (mas s 
density po oc r~ 2 ) or interstel lar medium (po ~ constant): see Ivan Paradijs et all (|2000h : 



Piranl (|2005l) : lMeszarod (|200d) for reviews. Light curves fluctuate strongly at early times, 



probably because of unsteadiness in the source; later, brightness falls approximately as 
a power law in time but often with undulations that may be due to inhomogeneities 
ahead of the shock. The observed radiation appears to be synchrotron emission, which 
implies that > 1CP 2 of the postshock energy density takes the form of magnetic field and 
highly relativistic electrons. Even after compression by the shock, typical circumstellar 
or interstellar fields would be many orders of magnitude too small. Therefore, it is often 
suppo sed that magnetic energy is created rapidly by plasma instabilities at the shock 



front (jMedvedev and I. oeUI !!)!)!)■. We wish to explore whether the compressed preshock 



field might instead be amplified gradually by macroscopic fluid turbulence. The source 
of the turbulence is supposed to be vorticity produced as the shock passes over inho- 
mogeneities in the ambient medium. Astrophysical applications, however, are deferred 
to a later paper. Our purpose here is to develop and test suitable theoretical tools: a 
relativistic version of geometric shock dynamics (hereafter GSD); and, independently, a 
suitable redefinition of relativistic vorticity that leads to conservation of the circulation 
on any fluid contour that does not cross a shock. 



The elements of nonrelativistic GSD were developed in the 1950s (MoeckeJ 11952 
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Chester! Il95i IChisnel]|ll957t IWhithair]|l957L Il958l Il959l ). IWhithaml <|l974h gives a ped- 



agogical review, upon which we have relied heavily. The one-dimensional version of the 
theory gives a functional or even algebraic relationship (rather than a partial differential 
equation) between variations in the pre-shock density and variations in the shock Mach 
number. The multidimensional version describes the effect of changes in shock area — 
divergence or convergence of the shock normals — on the Mach number. Thus, the theory 
reduces the dimensionality of the problem by one: in three dimensions, for example, it 
gives a closed set of equations for the evoluti on of the shock surface . GSD has even been 
adapted to reacting flows (detonation waves: iLi and Ben-dor 1 19981) . 

Naturally, there is a price to b e paid in accuracy for these simplifications. For a recent 
critique, see iBaskar and Prasad! (|2005l ). Nevertheless, GSD often performs remarkably 
well when there is reason to expect that fluid gradients or geometrical constraints near 
the shock ought to dominate, rather than reflections from boundaries behind the shock, 
and even in some cases where there is no such expectation. GSD successfully describes 
diffraction of shocks around corners, acceleration of converging shocks, and even the 
propagation of kinks ( "shock shocks" ) along shoc k fronts, as .judged by comparisons with 
experiment and with exact self-similar solutions (jBrvson and Grosslll96ll : ISchwendemanl 
1988; IWhithamlll974 and references therein). 

Rclativis tic unit s in w hich the speed of light c = 1 will be used. We adopt the con- 
ventions of Schutz ( 1990| ) for tensors; in particular, the metric in Minkowski coordinates 



= (x° = t, 



,x ) is rfpv = diag(-l, 1,1,1), while T^ v and T^ v denote the compo- 



nents of the same tensor in two Lorentz frames O, O. In all cases considered here, the 
energy density of the pre-shock fluid will dominated by rest mass, so that pressure and 
turbulent motions can be neglected ahead of the shock. 



1. Planar shocks 

Our goal is to transcribe GSD for an ultrarelativistic ideal fluid. Following Whitham j 



(1974), we begin with the case that the area of the shock is constant and the pre-shock 



density (po) is stratified on planes parallel to the shock front. In place of Mach number, 
we will be concerned with shock Lorentz factor (L) or rapidity parameter ($), the two 
being related by T = cosh3> ss e*/2 3> 1. These quantities are defined in the rest-frame 
of the pre-shock medium. 

The construction of GSD proceeds in two parts. First, the jump conditions are derived 
from the basic conservation laws; these relate the post-shock fluid properties to the pre- 
shock ones if T is given. This step is potentially exact but simplifies after approximations 
based on L 3> 1, p ~ 3P, and Po <C po- Next, characteristic equations are derived for the 
post-shock flow, and the (uncontrolled) approximation is made that one of the Riemann 
invariants has a known and uniform value behind the shock. 

1.1. Jump conditions 

For a planar shock propagating in the x 1 direction, the relevant components of the 
energy-momentum tensor are 

T 00 = (p + P) 7 2 - P, T 01 = (p + P) 7 2 f3, T 11 = { P + P) 1 2 I3 2 +P, (1.1) 

where the fluid 4-velocity has components 

Iff* (7, 7/3, 0,0)= (cosh <j), sinh 0, 0, 0) 

measured in the shock rest frame. The proper energy density p and pressure P are defined 
in the local fluid rest frame, so that they are Lorentz invariants. 
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The jump conditions in the shock rest frame are that T 01 and T 11 should be continuous. 
Since P as p/3 behind the shock and Pq <C pa in front of it, these conditions become 

T 01 : 4P 7 2 /? = -p Q T 2 T 11 : 4P 7 2 /3 2 + P = p Q T 2 . (1.2) 

Consistent with the ultra-relativistic approximation, the pre-shock three-velocity has 
been set to —1, which incurs an error ~ 0(T~ 2 ). Eliminating T 2 po between the two 
equations, dividing through by P, and setting 7 ~ 2 — > 1 — /3 2 yields (3/3 + l)(/3 + 1) = 0. 
The root (3 = — 1 corresponds to no shock at all. Therefore, (3 = —1/3 in the shock frame. 
In terms of the rapidity parameters of the fluid and the shock, tanh($ — 0) = 1/3. The 
latter is a covariant formulation since a Lorentz boost along x 1 with velocity v simply 
adds tanh -1 v to both and 4>. Substituting (3 = — 1/3 and 7 2 = 9/8 into either of 
eqs. (|1.2[) yields P = (2/3)r 2 p - So, the jump conditions are 

= $-tanh -1 ± = $-ln v / 2, (1.3a) 

C= f lnFK ^ (lnp + 2$-ln6) . (1.36) 

The peculiar factor y/3/4 will simplify the characteristic equations below. 

1.2. Whitham's Characteristic Rule 

Equations (| 1 . 3[) give two relations among the four variables (0, $, £, In po); or equivalently, 
(7, r, P, po)- They are exact up to terms 0(T~ 2 ). Whitham's formulation of GSD adds 
one more condition: the Riemann invariant associated with the characteristics that go 
upstream from the postshock flow toward the shock, P + , is supposed to have the same 
value as it would if the shock were a transition between constant states. The rationale is 
that the perturbations to the shock front are supposed to be localized; they are driven by 
small-scale density variations in the pre-shock fluid, or by local geometrical constraints 
on the shock, and therefore these perturbations are supposed to average out downstream. 

So the next step is to derive the Riemann characteristics fro m the equations of motio n 
T^ v u — 0. This has been done in greater generality elsewhere ( Marti and Miillerl 1 1 9 94h . 



but for completeness we shall rederive the special case we need. With eqs. (jl.ip . the 
equations of motion become 

T° v v oc (2 cosh 20+ l)-pC+ (4sinh20)</>+ (2 sinh 20)^=C' + (4 cosh 20)0' = 0, 
V 3 v 3 

T lv u oc (2 sinh 20) (4 cosh 20)0+ (2 cosh 20- l)-pC' + (4sinh20)0' = 0, 
V 3 v 3 

in which the dots denote 3/3x and the primes 9/Qx 1 . Rather than manipulate these 
equations directly, it is easier to boost into the local rest frame where = 0, find the 
characteristics there, and then boost back. Since cosh 20 — ► 1 and sinh20 — > 0, the 
equations above reduce to 

By adding and subtracting these, one sees that the characteristic velocities are ±l/\/3, 
and the corresponding invariants £ ± 0. Boosting along x 1 to any other frame simply 
adds a constant to 0. Therefore, 

R + = C, + is constant on C + : (^jf ) + = tanh ^0 + tanh -1 , (1.4a) 
P_ = ( - is constant on C_ : (|f)_ = tanh ^0 - tanh -1 ^=J . (1.46) 
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In Whitham's approximation, i?+ is constant not only along the C+ characteristics 
but everywhere in the postshock flow, even immediately behind the shock. Therefore, 
evaluating <f> and £ from the jump conditions (|1.3p . one obtains 

3 

$ + A lnpo ~ constant, where A = V3 - - » 0.232 (1.5) 

This approximate equation predicts how the shock speed slows in response to a transitory 
increase in pre-shock density: T oc p^ x . With eqs. (|1.3p . we get the corresponding changes 
in postshock rapidity and pressure: 

cf) + A lnpo = constant; P oc pJ _2A - (1.6) 

1.3. Comparison with an exact self-similar solution 

To reiterate, the approximation (jl.5p is intended to describe localized and transitory 
fluctuations in the propagation of a shock that has some prescribed average Lorentz 
factor r and advances into a "cold" medium with some prescribe d, but spatially va riable, 



pre-shock density p and negligible internal motions and pressure. IWhithaml (|1974T ) shows 
that the original nonrelativistic version of his theory approximates rather well the self- 
similar propagation of a planar shock from x < into a power-law density profile p oc 
(— x) n > even though this situation does not entirely satisfy the ass umpt i ons o f GSD. 
The corresponding ultra-relativistic solution has been given by ISari (|2006h . In Sari's 



terminology, the case of interest is a planar (dimensionality parameter a — 0) "Type II" 
shock with density exponent k = —n < 0. The Type I (II) solutions are those in which 
the powerlaw scaling of shock position with time can (cannot) be deduced from global 
energy conservation. Type II, where the scalings are determined by local conditions near 
the shock — rather than the inertia of the "piston" behind it — is the case for which one 
might hope that Whitham's theory would have some success. Indeed, Sari's equation 
(26) shows that the shock evolves as T oc p~ x with A exactly as in eq. (|1.5|l ! It is not 
clear why the agreement should be exact, but presumably the extreme simplicity of the 
ultrarelativistic fluid equations is somehow responsible. 



2. Non-planar shocks 

No vorticity can be created by an exactly planar shock, yet the ID theory above may be 
adequate for estimating the vorticity produced by encounters between an ultrarelativistic 
shock and a density inhomogeneity. Lorentz contraction causes the lateral dimension of 
inhomogeneities viewed in the shock or postshock frame to be larger by a factor r>l 
than the longitudinal ones, so that changes in speed and pressure are impressed upon the 
immediately postshock flow before it "notices" that the changes differ at other lateral 
positions. Thus it should usually be sufficient to evaluate the flow changes from the ID 
theory, and then take lateral derivatives to evaluate the resulting vorticity. 

Nevertheless, it is worthwhile to extend the ultrarelativistic version of Whitham's 
theory to nonplanar shocks for several reasons: 

• in order to study the stability of the shock; 

• in order to compare with exact spherical and cylindrical self-similar solutions, and 
with numerical tests such as refraction around an (oblique) corner; 

• because the extension is not difficult. 

The idea of Whitham's nonplanar extension is to insert a factor representing changes in 
shock area into the conservative form of the fluid equations. Thus let x 1 be a coordinate 
measuring arc length along the shock normal, and x 2 and a; 3 be coordinates in the shock 
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surface defined in such a way that a point moving along the shock normal maintains 
constant (x 2 ,x 3 ). The equations of motion are taken to be 

T^o + A- 1 (AT 01 ) 1 = 

T 10 + /T 1 (AT 11 ) x =0. (2.1) 

Here A is the 2D Jacobian relating the area of an element of the shock surface to its initial 
area. More precisely, since the characteristic equation does not hold across the shock, A 
represents the cross-sectional area of a bundle of streamlines immediately behind the 
shock; since the pre-shock medium is assumed to be at rest, the flow behind the shock is 
normal to it. 

Equations (|2.ip are not fully equivalent to T^ v . v = 0: for /i = 1, they do not contain 
the part of the covariant derivative associated with turning of the shock normal. They 
do however represent the divergence or convergence of the normals, which leads to area 
change and strengthens or weakens the shock. Because of the terms involving A, the 
quantities R± are no longer invariant along their respective characteristics C±. The 
equation for R + works out to 

d \ /Jm , . sinh 6 ( d \ , . 

— (( + &) = = — In A 

as J + V3 sinh + cosh 4> \ as J + 



' 'A] bxA, -2.2; 



+ 



where (d/ds) + = 3j + (v+) di is the derivative along the C + characteristic, and v + = 
tanh[</) + tanh _1 (l/\/3)] is the characteristic velocity. The final form of eq. (|2.2|) is in the 
pre-shock rest frame where cosh0 w sinh</> w 1 up to 0(T~ 2 ). 

The jump conditions (|1.3p are unchanged. Inserting these into (|2.2[) yields 

-i($ + Alnp + Mln^) »0, (2.3a) 
ds 

A=\/3-^, ^ = 3%/3-5. (2.36) 

Following Whitham, we have made the approximation that the characteristic equation 
applies on the shock, although its propagation speed [rapidity $ = </)+tanh _1 (l/3)] is not 
quite the same as that of the characteristic [tanh -1 v + — 4> + tanh - (1/V3)] Consistent 
with this approximation, the derivative d/ds in (|2.3p is taken to be the derivative with 
respect to arc length along the shock normal. Equation (|2.3[) predicts that the shock 
decelerates locally where its area increases, and accelerates where the area decreases. 
This will tend to stabilize corrugations in the shock front. 

Equations (|2.3p need to be supplemented by a vector equation for the shock normals. 
Introduce a function t(x, y, z) such that the locus of the shock in Minkowski coordinates 
is described by t = r(x,y,z). The normal to the shock is then n = Vr/|Vr|, and its 
3-velocity is V — n/|Vr|. The area function A of the shock satisfies 

V.(|)=0. (2.4) 



This is a purely geometrical, rather than dynamical, statement. I Whithaml (|1974 ) demon- 
strates it by applying Gauss's Law to a "flux tube" whose sides are made up of integral 
curves of H, and whose ends are elements of the shock surface at different times. 
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2.1. Comparison with nonplanar self-similar solutions 



Saril (|2006f )'s equation (26) for Type II solutions is equivalent to 



Alnr = a(5-3V3) - (V3 -f] ^-Inp, 
or \ 2 J dr 

p cx r _,c being the pre-shock density. Here a = 0, 1, 2 for planar, cylindrical, and spher- 
ical shocks, respectively. Since the area factor should scale as A cx r a in these three 
geometries, equation (|2.3|) predicts Sari's result perfectly. 

In these self-similar solutions, d imens ionality plays a limited role since individual shock 
normals are constant. Whit hair] (1974) discusses applications of the nonrelativistic the- 



ory to non-self-similar and truly multidimensional problems such as refraction of shocks 
around corners and obstacles. Here we can expect eqs. (|2.3l fe (|2.4|) not to be exact since, 
as noted above, they do not incorporate the full covariant derivatives in the equations of 
motion, and since the streamlines behind the shock are not perfectly straight. 



2.2. Detailed treatment of initially planar shocks in two dimensions. 

These details will facilitate comparison with numerical solutions of the full hydrodynamic 
equations for two-dimensional test cases (S}4|). They also serve to illustrate effects that 
involve changes in the shock normal, including the transverse propagation of disturbances 
along the shock front. 

We take z to be the ignorable coordinate. Following Whitham again, let ip be the 
angle between the normal and the x axis, so that dr/dx = \Vr\cosip = V^ 1 cos ip and 
dr/dy = V^ 1 sim/j. Equation ()2.4|) can then be rephrased as the two first-order equations 



3 3 

— (V 1 sinip) - — (V 1 cosip) = , (2.5a) 
ox oy 

3 3 

— (A- 1 costp) + — (A- 1 sin?/;) = , (2.56) 
ox v ' oy ' 

of which the first is simply the statement that d 2 r/dydx = d 2 r/dxdy. Together with 
(|2.3p . equations (|2.5p form a hyperbolic system. This is especially clear in the paraxial 
approximation where ip ~ 0(T~ 2 ). To this order, we may then replace simp — > ip, 
cos-0 — > 1, and V — > 1 — (2r 2 ) -1 , so that eqs. (|2.5[) become 

JU-. + 1. (_4-V) = 0. (2.61,) 

With our ordering, the term in ip 2 is of higher order than the others — it results from 
taking cos ip = 1 — ip 2 /2 rather than unity in the first of eqs. I|2.6p — but it does no harm 
and in fact makes the characteristic velocities work out more neatly. 

If the shock is initially planar and the pre-shock density initially uniform, then (j2.3[) 
implies that <!> + A In po + /j, In A is constant throughout the flow. Thus T in the first of 
equations (|2.6|) is to be regarded as a function of A and [x,y), given by 

r(A,x,y)=Tp x A»pv X A-», (2.7) 

in which the barred quantities are constants pertaining to the initially uniform medium 
and planar shock, and pa(x,y) is a prescribed function. The characteristic velocities of 
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the system (|2T6 j) - (j277]) are 

=^±^ 1/2 r- 1 . (2.8) 

y dxJ ± 

The factor of T _1 is easy to interpret as a consequence of relativistic beaming. A distur- 
bance propagating at the speed of light along an otherwise planar shock front would have 
transverse velocity dy/dt = ±T _1 in the pre-shock rest frame. To leading order in T -1 , 
we may replace dt with dx in this expression, so with /i 1 / 2 w 0.4429, the characteristics 
(|2.8[) are subluminal. 

Since eqs. (|2.5|) and (|2.6p are in conservation form, we may use them to study discon- 
tinuities in the shock front itself: "shock shocks." For a shock shock propagating at slope 
U = (dy/dx) ss , the jump conditions implied by (|2.6[) are 

[-2Utp + ip 2 - T- 2 ] = 0, (2.9a) 
[A-\^-U)]=0, (2.96) 

where [Q] denotes the discontinuity in quantity Q across the shock shock. Let us assume 
a homogenous pre-shock medium, po = p =constant, T ^constant. Then if tp = and 
A = 1 ahead of the shock shock, the post-shock-shock quantities (tp',A') satisfy 

(C/r)2 = W^T' V = U{1-A'). (2.10) 

Thus in the limit A 1 — > 0, we have tp' = U = 1/T, and it follows from eq. (|2.7p that 
T' — > oo. In the opposite limit A' > 1— but still A' < so that T' > 1 (else the 
ultrarelativisitic approximation would not apply) — we have %/)' = —{A')^/T = —1/T' and 
U = -ip'/V. 

Finally, because gamma-ray-burst shocks are believed to emanate from effectively 
pointlike explosions, it is of interest to consider a nearly spherical rather than planar 
shock. This case is effectively two-dimensional if the perturbations are axisymmetric. We 
take polar coordinates (r, 9, <j)) such that 9 = 0, 7t is the axis of symmetry and define tp 
to be the angle between the normal and radial directions, i.e. n ■ f — cos?/'. The analogs 
of eqs. (12. 5|) then become 



3 / r sin ip \ 3 / cos tp 



3r V V ) 30 V V 

3 f r 2 costp\ r 3 fsiri9sintp 



dr V A J sin0 30 V A 
and for tp 1, r 3> 1, eqs. (|2.6[) become 



= 0, (2.11a) 
= 0, (2.116) 



^WO + ^ 2 -r- 2 )=0, (2.12a) 



3r w ' 2 30 
3/r 2 \ r 3 / tp sin 



dr \ A J sin0 30 V A 

while (|2.7|) is unchanged. 



(2.126) 



3. Relativistic vorticity 

This discussion in this section is independent of the approximations of GSD, although 
the ultrarelativistic equation of state P — p/3 figures prominently. 
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As explained above, we are motivated by the need to explain the amplification of 
magnetic field behind the shocks associated with gamma-ray bursts, and by the possible 
role of turbulence in this amplification. Therefore, it may be worthwhile to record our 
assumptions about the relation of post-shock vorticity to the magnetic field. 

It follows from the induction equation of ideal magnetohydrodynamics, 

d t B = V x (vx B), (3.1) 
that magnetic energy increases according to 

J B-Bd?x = J BiBjdiVj d 3 £ , 

which involves the instantaneous shear (diVj + djVi — |V • v) and convergence (V • v) 
of the velocity field rather than the vorticity, which is its curl. Nevertheless, vorticity 
is important to secular amplification of the field by localized disturbances. In an ideal 
fluid, a localized nonvortical disturbance evolves into sound waves, whose oscillations 
produce only transitory changes in magnetic energy, and which propagate away from their 
source. Energy in vortical motions, however, remains localized, and the shear between 
neighboring eddies is expected to amplify the field exponentially on their turnover time. 



3.1. Vorticity and circulation 

Non-relativistically, the vorticity lu = V x v, where v is the fluid three-velocity. In a 
compressible but isentropic fluid without shocks, Kelvin's Circulation theorem is 

— <£v-dl = 0, (3.2) 

c 

where C is closed contour advected by the flow. 

The generalization of Q and eq. (|3.2|) to relativistic flow is not entirely straightfor- 
ward (|Eshraghj|2003[h Let U = (Z7° , J7 1 ,U 2 ,U 3 ) be the 4-velocity of the fluid, so that 



n^U^U" = — 1. In terms of the local rest-frame energy density p and pressure P, the 
energy-momentum tensor is 

T^ v = (p+ P)U»U V + rf v P. (3.3) 

The equations of motion 

T M "„ = (3.4) 

must be supplemented by an equation of state. Normally this involves two independent 
thermodynamic variables, e.g. P = P(p,T), P(p,N), or P(N,S), where T is the rest- 
frame temperature, N is the proper number density of conserved particles, and S is the 
entropy per particle. An essential feature of ultra-relativistic shocks is that the postshock 
particles are highly relativistic in the fluid rest frame, so that P — p/30 To the extent 
that the fluid is ideal [eq. (|3 - 3[) ] . entropy and temperature gradients then have no effect 
on the flow — except at shocks, but even there they do not have to be addressed explicitly. 

To illustrate this point, we consider a general equation of state in which entropy does 
influence the dynamics. The conservation of particle number is expressed by 

(NU^)^ = 0. (3.5) 



f This assumes that the stress is isotropic, which is not at all obvious in astrophysical appli- 
cations where the plasma is collisionless. 
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The vorticity turns out to be best formulated as 

^/nw = —Hp,v + Hi>,n > (3-6) 
in terms of the relativistic enthalpy h and its associated current H^: 

h = , H" = hU" . (3.7) 

With the First Law d(p/N) = TdS-Pd(l/N) in the form rT x dP = dh-TdS, eq. ([3^4]) 
can be rewritten as 

U u H^=TS^-h^. (3.8) 
The vorticity equation then follows from the "curl" of this, namely 

= T a S,f 3 -TpS, a . (3.9) 

The last line above vanishes if the entropy is uniform, S tll = 0, or more generally if there 
is only one independent thermodynamic quantity so that S = S(T). The intermediate 
expression is the Lie derivative of the vorticity considered as a 2-form, 

ft = dH = H v udx^ A dx v , 



with respect to the 4-velocity U. The statement of conservation of circulation is then 

d 

d7 



H a dx a = C G 1 1 fl = , 



dA A 

where A is a surface advected by U and dA is the curve bounding it, and d/dr = U^d^ 
is the convective derivative. Since different parts of the fluid move with different Lorentz 
factors, the surface A and curve dA will not (in general) remain within hyperplanes of 
constant Minkowski time t, unfortunately. 

Now we specialize to the ultra-relativistic equation of state P = p/3. The true entropy 
per particle is 5 oc ln(P/iV 4 / 3 ), which will not be uniform after the shock passes over 
density inhomogeneities. However, if we define an ersatz number density N oc P 3 / 4 
at some initial time, so that P/7V 4 / 3 is spatially uniform, and if we demand that N 
evolve according to (|3.5|) with N instead of N, then P/N 4 ^ 3 will remain uniform in 
smooth parts of the flow, though not across shocks. This follows because Uf/T^" v = 
-W v P lV - iPU^ = 0, whence U V {P/N 4 / 3 )^ = 0. 

We don't actually have to deal with TV directly. Since h = (p + P)/N = 4P 1 / 4 , we can 
simply redefine the enthalpy current as 

H n = p i/i u n (when P = p/3 only). (3.10) 

Then eq. f|3.8|) becomes 

C/^,, = -(V /4 )^ P=(H»H fl ) 2 . (3.11) 

The vorticity defined in terms of this H via (|3.6p is conserved in the sense that the 
righthand side of Q3.9|) vanishes in all smooth parts of the flow, even after shocks. In 
particular, if £1 = initially then it remains zero as long as the flow remains smooth. 
The jump conditions that follow from integrating (|3.4[) across shocks are not equivalent 
to the corresponding integral of (|3.9p . however, and so this vorticity can be created at 
shocks. 
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3.2. Vorticity vector 

fl a p — —flpa would appear to have six algebraically independent components. In fact it 
has only three. To see this, note that 

-h,„ = {H V U V )^ = H Vtll U v + H V U V ^ = H V ^XJ V + h{U v U v ) iP = H V ,^U V . 

Therefore if the fluid is isentropic, or if it is ultrarelativistic and the enthalpy current is 
defined by eq. (|3.10p , then eq. (|3.8[) becomes 

U v n^ = 0. (3.12) 

This means that in the local rest frame, where U v — ► Sq, the "electric" components 
Q i0 = — r2 i of the vorticity vanish, and only the three "magnetic" components flij = 
—Fiji survive. Let Co be the three-vector field with components u>i — eijk^ljk/^ m an 
arbitrary inertial frame (not necessarily coinciding with the local fluid rest frame) . If one 
uses eq. (|3.12[) to eliminate the inertial components £lio from the identity 

which is the tensorial expression of d(diJ) = 0, the result is 

d t uj - V x (v x w) = 0, (3.13) 

where v is the three- velocity, v l = U l /U . This is formally identical to the nonrelativistic 
vorticity equation of an isentropic fluid, except that Co is V x H rather than V x v. The 
derivation just given, which parallels that of the induction equation (|3.ip . shows that 
eq. (|3 . 13(1 is relativistically covariant. 



4. Numerical simulations and tests 

In order to test the predictions of the GSD theory we have performed a series of numeri- 

cal sim ulations with the the RAM special relativistic hydrodynamics code lZhang and MacFadven 
( 20061) . The simulations test the variations of pressure and Lorentz factor behind a strong 
shock as it passes over density perturbations in the pre-shock medium. The simulations 
are performed on the domain x = [0.0,1.0] with reflecting boundary at x = 0.0 and 
outflow (zero gradient) boundary at x = 1 and adiabatic index j a = 4/3. In the notation 
of the present paper, the equation of state is P = — l)(p — mN), where m is the rest 
mass per particle. 

We have performed ten simulations in which the pressure is initially set to P = 10 5 
for x < 0.001 and P = 10~ 6 elsewhere, to create a plane-parallel analog of an ex- 
plosively driven spherical gamma-ray-burst shock. The rest-mass contribution mN to 
the energy density for x = [0.8, 0.9] is initially set to one of the ten perturbed val- 
ues p p = {1.0, 1.2, 1.4, 1.6, 1.8, 2.0, 5.0, 10.0, 20.0, 30.0} and to Pl = 1.0 elsewhere. The ten 
simulations are otherwise identical. A strong relativistic shock with Lorentz factor r ~ 10 
is initially driven into the medium as the over-pressured region expands and a thin rela- 
tivistic shell forms behind it. As the shock crosses the (positive) density perturbation, its 
Lorentz factor decreases and post-shock pressure increases with the values seen in Fig. 
1. Note that 7 is the Lorentz factor of the fluid rather than that of the shock itself; for 
an ultrarelativistic shock advancing into a cold, stationary medium, the jump condition 
(|1.3a[) implies that these two Lorentz factors are represented simply by 7 = a 
relation that does not depend upon the approximations of GSD. We measure 7 and P 
at t = 0.9, after the shock has crossed the perturbation region and is at x « 0.9, and we 
record for that time slice the values P m and 7 m where the pressure reaches its maximum 
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Figure 1. Comparison between analytic GSD predictions and numerical simulations for 
post-shock Lorentz factor and pressure as a function of the pre-shock density perturbation. 
The right panel shows the fractional difference between the GSD prediction Xgsd and the sim- 
ulation value X s i m , where X stands for maximum pressure P m (plus signs) or Lorentz factor 7 m 
(asterisks). 



post-shock value. In Fig. 1 (left panel) we show the logarithms of P m (plus signs) and 
7~ (asterisks) versus \og(p p / Pi). The ordinates have been scaled by the values P and 
7o obtained from a fiducial run with uniform density (p p — pi — 1). The lines show the 
scaling predicted by the GSD approximation with slopes A = \/3 — § « 0.232 (solid) 
and 1 — 2 A ~ 0.536 (dashed). In Fig. 1 (right panel), we show the fractional difference 
between the GSD prediction and the numerical simulation for maximum pressure (plus 
signs) and Lorentz factor (asterisks). 

We find that the GSD approximation works surprisingly well even for large density 
perturbations, justifying its use up to density contrasts of factors of ten or more. 

We acknowledge the use of the Scheides cluster at the Institute for Advanced Study. 
The software used in this work was in part developed by the DOE-supported ASCI/ Alliance 
Center for Astrophysical Thermonuclear Flashes at the University of Chicago. 
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